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Abstract 



We study the model of (2 + l)-dimensional relativistic fermions in a ran- 
dom non-Abelian gauge potential at criticality. The exact solution shows that 
the operator expansion contains a conserved current - a generator of a contin- 
uous symmetry. The presence of this operator changes the operator product 
expansion and gives rise to logarithmic contributions to the correlation func- 
tions at the critical point. We calculate the distribution function of the local 
density of states in this model and find that it follows the famous log-normal 
law. 
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Typeset using REVT^ 



I. INTRODUCTION 



In this paper we continue to study the model of (2 + l)-dimensional massless Dirac 
fermions interacting with a random static non-Abehan gauge potentiaL The Hamiltonian 
(or rather a generating functional for the Green's functions) of this model is given by 



iH-eJ = J d'x {Rt,Lt 



\ -en (5x + i^J - y) ) 

where the random fields ^^^{x^y) (a = ±) have a Gaussian distribution: 

(A^^(fi)A^^(f2)) = A5{n - f2)5,,„b(5,^5^^ (2) 

The fermionic fields Ra-, La represent respectively the right and the left moving components 
of the spinor field, and a takes the values 1, A^. This model was introduced in Ref. to 
describe normal excitations in two dimensional non-s-wave superconductors with disorder. 
In this context denotes the number of nodes of the order parameter on the Fermi surface. 

Since the disorder is time-independent we consider the Fourier components of the 
fermionic fields with different frequencies separately. This reduces the dimensionality of 
the problem, making it two dimensional. The model (1) is exactly solvable in the subspace 
e = 0, where it is critical. The presence of the superconducting order parameter fixes the 
chemical potential and thus insures that the disorder is diagonal in chirality. This is es- 
sential for the criticality at e = 0. At the critical point, one can apply the methods of 
conformal field theory and calculate scaling dimensions of the fields and of their multi-point 
correlation functions. This gives us a rare opportunity to obtain nonperturbative results for 
a non-trivial random theory. We hope that a study of this exact solution will give an insight 
into general properties of random systems. 

The averaging over the disorder can be done either through the replica trick or using 
the supersymmetric approach , . In this paper we shall mostly use the replica approach 
with which we are more familiar. To demonstrate that the two approaches are equivalent we 
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shall (i) compare the conformal dimensions of the primary fields and (ii) demonstrate that 
both representations give the same conformal blocks for the four-point correlation function. 

As in the standard localization theory (see for example 0), one can integrate out the 
fast degrees of freedom and derive an effective action for the slow ones in the form of a 
sigma model. This program was carried out in Ref . [|I| . The resulting sigma model has the 
following action: 

S = So + MenJ d2xTr(Q + Q+) (3) 

where the Sq action contains the Wess-Zumino term: 

Nr 



\ T C 

So= -Y J <^Md,^Y + NW[SUir); g] 



(4) 



Iff 2 

W[SU{2ry,g]= — J d\ \TT{d,g+d^g) + - d^eabcTrig+dagg+d^gg+d^g) 

where Q is the r x r matrix Q = g exp[i\/47r$], and where g belongs to the SU(r) group. $ is 
a real scalar field defined on the cirle with circumference a/tt. The quantity M is the energy 
scale introduced by the disorder: M ~ exp[— 27r/A^^]; it marks the crossover from the bare 
density of states (DOS) p(e) ~ |e| to the renormalized DOS p(e) ~ {el" {u = {2N^ - 1)-^). 
M serves as the ultraviolet cut-off for the sigma model (^. 

It is well known that in conventional localization theories the average DOS is not affected 
by the disorder (its higher moments are, however). It is not the case for the model (1), where 
the DOS is directly proportional to the order parameter. As it was shown in |T|], the local 
DOS is given by 

p{e,x) = ^TT[Q{x) + Q+{x)] (5) 

This means that even the average DOS is strongly renormalized. At e = the sigma model 
(^) is critical and (p(0, x)) = 0, as it might be expected in a critical theory in two dimensions 
0. In three dimensions a finite DOS emerges at e = ||^. In our previous publications we 
interpreted this effect as a manifestation of violation of some continuous symmetry present 
in the theory ( |T|, 0). However, the meaning of this symmetry has remained obscure. In 
this paper we identify the operators which generate this symmetry and derive their algebra. 



The paper is organized as follows. In Section II we discuss general properties of the model 
at criticality and derive the expression for the conformal dimensions of its primary fields. 
In Section III we derive the expression for the four-point correlation function of local DOS. 
It turns out that this correlation function contains logarithmic singularities and therefore a 
fusion of local DOS generates operators with unusual properties - the so-called logarithmic 
operators. In Section IV we develop a general theory of such operators and demonstrate that 
the appearance of such operators always implies the presence of some continuous symmetry. 
It may well be that this symmetry is present in all critical models with disorder. In Section V 
we show how these general results hold for the model (4). In Section VI we study correlation 
functions in the model deformed by the logarithmic operators away from criticality. It turns 
out that even in the case of a marginal deformation the correlation functions have logarithmic 
corrections. In Section VII we derive the log-normal distribution function of the local DOS. 
The paper contains a conclusion and an appendix where the conformal blocks of replica and 
supersymmetric theories are compared. 



II. CONFORMAL DIMENSIONS 

The WZNW model has been well studied at finite r. There is an extensive literature 
on the subject, but we particularly recommend the original publication by Knizhnik and 
Zamolodchikov |]T2[. These authors derived explicit expressions for the four-point correlation 



functions of primary fields which we are going to exploit. In our calculational procedure we 
follow the general principle: when calculating any ra-point correlation function ^^(1,2, ...n) 
r is treated as an arbitrary number on all intermediate steps of the calculations until the 
final expression is obtained. We define the replica limit as follows: 

2N 

F(l,2, ...n) = lim F^(l,2,...n) (6) 

The reason for the introduction of the extra factor 2N will be discussed later. 

Let us study the correlation functions of the Q, Q^-fields. The problem of indices is 
simplified by the fact that the Qpq matrices are slow parts of the operators 



N 

Qpq ~ RXpLa,q (7) 

From this fact one can derive a simple recipe for the index structure of n— point corre- 
lation functions: it is the same as for the ra-point function of the Yl,a -Ro p-^a.g-helds in the 
theory of massless free fermions. The simplest example is the 2-point function [0,0: 

where 1/2N'^ is the conformal dimension of the composite operator Q given by the sum of 
the dimensions of the bosonic exponent exp[iV^$] and of the operator field gpr from the 
fundamental representation of the SU(r) group: 

1 



A = lim 



1 (r — 1/r) 
2WV ^ 2{N + r) 



2m 



(9) 



In the replica limit we get from Eq.( 



G{zrz) ^ lim(2iV/r)(Tr[Q(^,^)]Tr[Q+(0,0)]) = {M\z\r"'' (10) 

All other operators are generated by fusion of the fundamental fields TrQ and TrQ^- The 
corresponding primary fields are composite fields of bosonic exponents and Wess-Zumino 
tensors belonging to irreducible representations of the SU(r) group. These representations 
are classified by Young tableaus which can be represented by a string of numbers /i > /2 > 
••• > /r > 0. Only representations with fi<N — l are generated |T^. The corresponding 
conformal dimensions are given by the expressions 



N + r 2rN 



Cf=lUf? + {r + l-2l)f,]-^ (11) 



where f = Y.i fi- the replica limit we get 



A/ = ^ + ^E[/N(2/-l)/z] (12) 
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This expression coincides with the conformal dimensions obtained by the supersymmet- 
ric approach after rows and columns in the Young tableau are interchanged. For in- 

m 

stance, for the representation (1, 1, ...1, 0, ..) we reproduce the expression obtained in P(see 
Eq.(4.48b)there) for the (m, 0, ...0) representation: 

(A^ - l)m 



m 
2N 



N 



(13) 



III. FOUR POINT CORRELATION FUNCTION OF THE ORDER PARAMETER 

FIELDS 

Let us now study the four point correlation function of the Q, fields. The index 



structure is the same for all r 12 



-1 



1^13^241 



(W + W), (14) 



\Zi2Zi4Z23Z34\ 

W= [^PW2^p-iPA^qiq2^qzqAWll{^^^) + (^PlPa '^P2P4 (^gigs '^9294 ^^22 (a^, x)] , (15) 

W= [5pip25p3P4^gig3^92g4^i2(a;, x) + 5p^p,^5p^p^5q^q^5q.^q^W2i{x, x)] (16) 



where 



^ _ ^12^34 _ _ ^122:34 ^yj^ 
Z13Z24 ^13^24 

Here the functions Wab{x,x), {A,B = 1,2) satisfy linear differential equations (the 
Knizhnik-Zamolodchikov equations) which we shall discuss later in detail. Now note that in 
our theory we shall deal only with correlation functions of TrQ, TtQ^ (for simplicity we do 
not consider transport phenomena which would need introduction of advanced and retarded 
correlation functions). Since all correlation functions must be proportional to r, only the 
W term (that is the term with all indices equal) survives in the replica limit, the W term 
being proportional to r^. Therefore we have 
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(2Ar/r)(TrQ(l)Trg+(2)TrQ(3)TrQ+(4)) 

-1 



2NM-*/^ 



1^13^241 



\Zl2Z34^Z23Z34,\ 



[Wi2ix,x) + W2iix,x)] (18) 



The functions Wab{x,x) = U^'^W^ {x)Wj^ (x) are composed of linearly independent 
solutions of the Knizhnik-Zamolodchikov equation, written W^\x) and Wl^\x) (conformal 
blocks). In the replica limit these equations have the following form: 

Nx^ = -W2, N{1 - x)^ = Wi (19) 

ax ax 

Thus for the function Wi we get the following hypergeometric equation: 

dx \ dx J 1 — X 

Here we encounter a problem. A hypergeometric equation has always two linearly indepen- 
dent solutions, normally expressed in terms of powers and hypergeometric functions. Usually 
there are three sets of solutions defined in the vicinity of x = 0, 1 and 00 respectively. These 
pairs of solutions are related to each other via simple transformation rules (see, for exam- 
ple [|T3|). Eq.(|20) is an exlusion: in the vicinity of x = one of the solutions contains a 



logarithmic singularity and cannot be expressed in terms of hypergeometric functions (this 
second solution was overlooked in the previous publication of one of the authors ) : 

Wi^\x)= F{l/N, -1/N, 1, x), 
Wi^\x)= \nxWi°\x) + Hi{x) 
NWi'^\x)= xF{l + 1/N, 1 - 1/N, 2, x), 

NWi^\x)= \nxW^^\x) -N^ + H2{x) (21) 
where Hi^2{x) are functions that are regular at x = 0: 
Hi{x) = 2^ [^(1/A^ + n)- ipil/N) + 

n=l 

+ij{-l/N + n) - ^{-1/N) - 2^{n + 1) + 2^(1)] 



CXD 



n=0 '^•v"' + 

+^{1 -1/N + n) -tpil- 1/N) - tpin + 1) + ^(1) - tp{n + 2) + ^(2)] (22) 



where (a)„ = r(a + n)/r{a). Only in the vicinity oi x — oo are the solutions still hypergeo- 
metric functions (for N ^ 2). 
At \x\ << 1 we have 

Wi°\x) = 1 + 0{x), Wi^\x) = \nx[l + 0{x)], 

NW^°\x) ^x + O(x^), NW^^\x) = -N^ + x\nx[l + 0(x)] (23) 

Now we have to choose the matrix C/^* in such a way that the resulting expression for the 
four point correlation function be a uniquely defined function in the complex plane of x. It 
also must be invariant under the permutation of points 1 and 3 (2 and 4) which means the 
invariance under (crossing symmetry) . These properties are achieved 

when 

um ^ ^(10)^ ^(11) ^ ^ ^^(01) (24) 

To find h, we first note that the solutions to the Knizhnik-Zamolodchikov equations obey 
the monodromy properties 

Wi''\l - x) ^ aiW^\x) 

Wi''\l - x) ^ biWi'\x) (25) 



where 



ao= ar[iP{l/N) + i^{-l/N) - i^{2) - V'(l)] 



r(i/A^)r(-i/A^) 

6o= ai[^{l/N) + i^i-l/N) - 2V'(1)] (26) 



Using the crossing symmetry, we then find that 



h = ^-^h±^ ^ 1/2 (27) 
aiOi 



Thus we get 
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G(l,2,3,4) 



2N 



\z13Z24\M ^ 

\Zi2ZmZ23Z34\ 



X [Wi''\x)wi^\x) + Wi^\x)wi''\x) + ^wP{x)wi''\x) + {x^ x)] (28) 



Here we choose U^'^^^ = —l/AN'^ for normalization. 

In order to derive the operator algebra of the model we consider various limits of this 
formula. In the limit ^43 = e ^ we get 



([gQ+(3)]g(i)Q+(2)) 



1 



£2:12 



12/ 



4f- 



(^Zi2 



+ c.c I In |e 



Z12 



(29) 



N'^ \Z13Z23 / V " Z13Z23 

From now on we shall use Q without subscripts instead of TiQ assuming that the replica 
limit has been taken. We shall also put M = 1. 

The three-point correlation function (p9D is very unusual from the conformal field the- 
ory point of view because it contains logarithms. Therefore we pause to consider general 
properties of logarithmic operators. 



IV. GENERAL PROPERTIES OF LOGARITHMIC OPERATORS 

So far correlation functions with logarithms at criticality have been obtained in the 
WZNW model on the supergroup GL{1, 1) [§, in the C = - 2 model ^ and in gravitationally 
dressed CFT ||T0|. It was first pointed out by Gurarie in Ref. that the appearence 



of logarithms in correlation functions is due to the presence of special operators, whose 
operator product expansions (OPE's) display logarithmic short-distance singularities. These 
logarithmic operators have conformal dimensions degenerate with those of the usual primary 
operators, and it is this degeneracy that is at the origin of the logarithms (cf. our discussion 
of the degenerate hypergeometric equation in the previous Section). As a result of this 
degeneracy one can no longer completely diagonalize the Virasoro operator Lq, and the new 
operators together with the standard ones form the basis of the Jordan-cell for Lq. In order 
to get a better insight in the situation, we shall consider the simplest example, which was 
mentioned in ||10i, namely, the Liouville model with the action 



(R^^^ is the Riemann curvature on a two-dimensional surface - the world sheet) with the 
stress-energy tensor 



and the central charge 



C = 1 + 3Q2 (31) 

The primary field exp («(/)) has a dimension 

A„ = a{Q - a)/2 (32) 

This means that there are two operators with the same dimension A^, namely V± = 
exp(a-|-0), where 

Q , 1 



If = 8Aq,, i.e. when a = Q/2, there is a degeneracy a+ = a_ and instead of two 
exponential primary fields we have only one exponent C = exp(|Q0) and the new operator 
D = (f)exp{^Q(f)) with the same dimension A = The latter field is sometimes called 

the puncture operator. It was discussed in in the context of the Liouville gravity, when 



the action (|30|) describes the gravitational (Liouville) sector of a (non) critical string in the 
conformal gauge. 

It is easy to get the OPE of the stress-energy tensor T with these fields. After simple 
calculations we find 

T{z)C{0)= 4c(0) + -d,C{0) + ... 

T{z)D{0)= ^D{0) + ^C(O) + ^d,D{0) + ... (33) 

where the dimension of the fields C and D is A = Q'^/8 and the normalization of the field 
D was defined as -D = (2/(5)0 exp(i( 
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It is easy to see indeed that there is a mixing between C and D and that the Virasoro 
operator Lq which is defined through the Laurent expansion T{z) = J2n LnZ~"'~'^ is not 
diagonal 

Lo|C>=A|C>, Lol/^ >= A|D > +|C > (34) 

Let us also note that usually one can think about factorization of the primary field into the 
product of chiral left and right operators using the decomposition (p{z,z) = (pL^z) + 0_r(^), 
leading to exp[a0] = exp[a0i^(2;)] x exp[a0/j(2)]. For a logarithmic operator we have 

(j) exp[a;0] = exp[a0L(z)] x exp[a0ij(z)] + (pR^z) &c^[a(j)L{z)] x exp[a0ij(^)] 

and thus the logarithmic operator is the sum of left and right operators each of which can 
be factorized. 

This simple example illustrates a quite general property of all theories with logarithmic 
operators and OPE (|33D is valid in all these theories. One can also obtain some general in- 
formation about two- and three-point correlation functions with operators C and D starting 
from the four-point correlation function ||^ 

{A{z,)B{z,)A{z,)B{z,)) = _ ^^^^^l^^^ _ [x{l - a;)]^-^-^-F(x) (35) 

where x is defined by Eg. ([T7|) , and where we have extracted the factor [x(l — a;)]^c;-AA-As 
to make F{x) finite at x — or x — > 1 in the ordinary case. In the case of logarithmic 
operators one will get F{x) = {d + clnx + o{x)) at small x. To reproduce the logarithmic 
singularity at x = after the fusion of A{zi) and -8(^2) one has to postulate the following 
OPE (we restrict ourselves to the chiral sector): 

A{zi) B{Z2) = (zi - ^2)^c-a.4-Ab + c\n{z, - z^) + ..] (36) 

Taking the limit Zi — ^ Z2 one immediately gets from the four-point correlation function the 
following three-point correlation functions: 

{C{zi)A{z^)B{zi))= Aa+Ac-As As+Ac-Aa Aa+a^-Ac 

^13 ^14 ^34 
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Now let us consider the ^(-23) and B{z4) fusion which after insertion of (|36D into 
will lead to the following two-point correlation functions: 

{C{x)D{y))= {C{y)D{x)) 



{x - y) 



2A 



c 



{D{x)D{y))= _ ^^^^^ (-2cln(x - y) + d) 

{C{x)C{y))= (38) 

The first equation imposes a strong constraint on the dimensions Ac of logarithmic 
operators, namely that the dimension Ac must be an integer. To prove it let us note that 
we have {C{x)D{y)) = {C{y)D{x)), which means that the correlation function is invariant 
under the permutation of x and y, which means that (— 1)^^^; = x go = n. In case of 
noninteger Ac one must have the structure constant c = and only the operator D will 
survive in OPE (|36D , however in this case it will be an ordinary, nonlogarithmic operator. 

This new result about dimensions of logarithmic operators means that for any logarith- 
mic operator we have a hidden continous symmetry. This symmetry is generated by the 
conserved holomorphic (or antiholomorphic) current C{z). This current is a symmetric ten- 
sor of rank Ac, which is a usual vector current if Ac = 1. In the next Section we shall 
demonstrate the existence of such a conserved vector current in the model with disorder that 
we are considering. Let us note that we have also proved that there is no central extension 
in the corresponding current algebra, or in other words there is no anomalous Schwinger 
term in the current- current commutator. This is the direct consequence of the triviality of 
the correlation function {C{z)C{0)) = 0. 

V. OPERATOR PRODUCT EXPANSIONS IN THE MODEL (4) 

In this Section we demonstrate how the general theory just discussed applies to the model 
(4). For this end we study OPE's in this theory after the replica limit has been taken. In 
doing so we assume that the replica limit can be described by some quantum field theory. 
Here we encounter a certain ambiguity, namely, that we can define the correlation functions 
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with an arbitrary pref actor. It turns out that this pref actor is fixed by the requirement of 
self-consistency of OPE's. The latter is achieved when one uses the definition (|^). This 
explains the necessity of the factor 2N in Eq.(^). 
We suggest the following OPE: 



Q{z)Q+{0) = 1^1^^/^' {l-z [d{0) + C(0) In \z\^] - z [d{0) + C{0) In \z\^] ...} (39) 



where D, C and Z), C are some new operators whose correlation functions are to be found. 
Notice that in the conventional WZNW theory this operator expansion would contain the 
unit operator and the operator in the adjoint representation. However, the latter one has 
the conformal dimensions which vanish in the replica limit: 

Therefore we have here a situation described in the previous Section: the conformal di- 
mensions of descendants of the unity operator become degenerate with the dimensions of 
descendants of some other primary field (the adjoint operator) which gives rise to logarithms. 
Substituting (|39D into Eq . (^Of) we get 

(g(l)Q+(2)C(3))=-i,|zi.r^/^^ 



(g(l)g+(2)D(3))=iV-2|;,i2r'/^'^^ln|^^| (41) 



2^13 2^23 2:13^23 



Setting Zi2 = e in these equations and using the OPE (|39|) we get the following set of 
two-point correlation functions: 

<^W^(2))^ 
(C(1)C(2))=0, 

(^(1)^(2))= (42) 

There are similar expressions for C*, -D-operators with z being substituted for z. Correlators 
of C, D and (7, D are equal to zero. 

Setting ^31 = e in Eqs.(0) we deduce the following OPE: 
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c^(^Mo)=^g(o) + ..., 

D{z, ^)Q(0)= In |z|Q(0) + ... (43) 

with the same equations for Q^, except for a change of sign. These OPE's and the fact that 
C{z) does not depend on z {d{C{z)D{0)) = 0), enable us to identify C,C as generators of a 
continuous symmetry. It is this symmetry which is associated with the order parameter p. 

Conformal field theories are characterized by their symmetry group and a number C 
called 'conformal charge'. Formally C is a coefficient in the pair correlation function of 
stress-energy tensor operators. A physical meaning of C becomes clear when we recall that 
a theory with an integer conformal charge C = k is equivalent to the theory with k species 
of free bosonic fields. Thus C in unitary theories counts an effective number of degrees of 
freedom. The central charge of our theory is the sum of central charges of the free bosonic 
field (C = 1) and the WZNW model on the SU(r) group: 

Thus the resulting central charge vanishes, as it must be; however, according to the definition 
of the replica limit (^) the physical correlation function of the stress-energy tensors remains 
finite: 

(T(.)T(0))=lim^^-^ = 4 (45) 

Superficially this looks like the effective central charge Ce// = 2. However, Cg// does not ap- 
pear in the fusion rules of the stress-energy tensor components inside of correlation functions 
with matter fields, where we have 

T(^)T(e) = + + - (46) 



As we have mentioned above, the numerical coefficient in ( ^Sj) is fixed by the self-consistency 
requirements of OPE. 

Applying twice the OPE ( pUD to the four-point correlation function and using the Ward 
identities for Q fields, we get the following set of identities: 
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{T{z)C{l)D{2))= Y: 



j=l,2 k ^j) 



Z — 



-d, (C(1)D(2)) 



(f(^)C(l)D(2))= -— a,(C(l)D(2)) = 



(T(z)z}(i)D(2))= 5: 

i=l,2 k V- -J 

1 



^ + {D{\)D{2)) + E ^(C(1)D(2)) 
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(T(z)Z}(l)Z}(2))= 



i=i,2 k ^ 



d, {Dil)Di2)) + 5: 



1 



j=l,2 ^iJ 



j=i,2 ^ 



■(C(1)D(2)) 



(47) 



We can then substitute the two-point correlation functions to get 
{T{z)C{l)D{2))- ^ ^ 



2iV2 (2 - - 22)2 

(T(z)Z}(l)D(2))= 4 ... ..J... .,,J lnki2|-l/4) 



(T(z)5(l)D(2)) = 



1 



"12 



(48) 



2iV2 {z - ziYiz - 22)2^12 
Taking into account Eq. (^5]) and Eqs. (^) we conclude that these expressions are compatible 
with the following OPE: 



2m 



+ T{0 + - 



D{z,z)D{^,0=j;^ 



\n\z-f\ (z-f)^ - - 

I ^1 +(ln|z-e|-l/4)T(0- ^L t{0 + 



Aiz-0' 



and 



[z - 0^ - _ - 



(49) 
(50) 

(51) 

(52) 

(53) 
(54) 



From the OPE's (^, ^9]), the Ward identity for the stress-energy tensor and primary fields 
and the Knizhnik-Zamolodchikov equation ( |T^) we derive the following Ward identity: 



{C{z,)D{z2,Z2)Q{^U^i)-Q^{^2N,^2n)) 



Er^(^(^2)g(ei,6)-Q+(6iv,6iv)) - ^^^(Q(ei,ei)-Q+(6iv,6iv)) (55) 



1 



2iV2 Y ^1 - 



15 



where a = 1 for Q and —1 for Notice that the operator D does not appear in the right 
hand side of this identity. This Ward identity is an important one since it, together with 
Eq. ( ^9]) estabhshes an isomorphism between the representations of the Virasoro algebra and 
the algebra of the conserved current C. 

Now let us study the fusion of Q with itself. For this end it is more convenient to rewrite 
the four-point correlation function (^) in terms of the solutions regular at x — * oo. For 
Eq.(|0D these solutions are (A^ ^ 2) 

Wi''\x)= (-x)-^/^F(l/Ar, 1/N, 1 + 2/N; 1/x), 
Wi^\x)= {-x)-^/^F{l + 1/N, 1 + 2/N; 1/x) 
Wl^\x)= {-xy/^F{-l/N, -l/N, 1 - 2/N; 1/x), 

Wi^\x)= -{-xf'^F{l - l/N, -l/N, 1 - 2/N; 1/x) (56) 
These solutions have extremely simple monodromy properties: 

-x)= w!f\x), w!f\i - x) = lyr(x) 

W[^\l - x)= -W^^\x), W!^^\l -x) = -W[^\x) (57) 
The crossing invariant form of the correlation function is 

W{x, x) = a[wl'^\x)wi'^\x) - k'^wl'^\x)W^^\x) + (x -> x)] (58) 

where 

r(l + 2/iV)r^(-l/Ar) 

r(i -2/iv)r2(i/Ar) 

The coefficient a whose numerical value we do not provide should be choosen to match 
Eq . (^8]) to the correlation function (^) regular at x = 0. Let us consider the limit 2:31 = 
e ^ we have 

G(l, 2; 1 + 0, 2 + 0) = 2a|e|-^/^'[|z/e|^(^-2)/^' - e\z/e\-'^''+'^/'''] + ... (59) 

This expansion is valid only for N ^ 2. In this case it corresponds to the standard operator 
product expansion: 
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g(l)Q(2) = Cy'|^i2r'^+'^^OA(2) + Cl/'|zi2|-'^+'^^05(2) + ... (60) 

where Ci, C2 are numerical coefficients and Oa and O5 are operators from the asymmet- 
ric and the symmetric representations whose Young tableaus are (1,1,0,...) and (2,0,...) 
respectively. Their conformal dimensions are given by Eq.(|l^): 

2-N ^ 2 + N 

which reproduces the result obtained in the previous publications [0] and [^. 

At = 2 the dimension of the antisymmetric operator vanishes. Now we have a situation 
where there are three operators with zero conformal dimension - the unity, the adjoint oper- 
ator and the operator in the antisymmetric representation. This situation will be discussed 
elsewhere. 

VI. DEFORMATION BY THE LOGARITHMIC OPERATORS 

The conventional WZNW model remains an integrable theory even if one changes the 



coefficient in front of the Tr(c}^(7"'"c}^(7)-term in the action (4). According to |T^, such per- 
turbation is equivalent to the J_i J-i'l'^^-operator (recall that is the primary field in the 
adjoint representation). The corresponding beta function is 

^(7) - ^7 (62) 

where 7 is the deviation of the coupling constant from its critical value. In our case c„ = r ^ 
and the beta function apparently vanishes. This means that the perturbation becomes 
marginal and we have to reconsider the terms of higher order in 7. Despite the fact that 
(^ab (-[Qgg T^Q^ appear now in OPE, its decendants, that is the logarithmic operators D 
do appear. We suggest that the change in the coupling constant of the WZNW model (4) is 
associated with the perturbation by the marginal operator 'jDD. We warn the reader not to 
confuse this perturbation with a change of the disorder strength A which is truely irrelevant, 
leading to a change of the cut-off M. One physical mechanism of a marginal deformation 
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away from criticality in the model (4) was described in [|l| (see Chapter 7). We conjecture 
that this deformation is generated by the 7Z)D-perturbation. 

In the case of a deformation 7 / d^zO{z,z) caused by a usual marginal operator O one 
has two possibilities depending on the operator product expansion 

O{z,z)O{0) = f^ + ... 

The first one is when / = 0, i.e. the OPE of 0{z, z)O{0) does not contain the operator 
O itself. In this case this operator is truly marginal and one has the continuous family of 
conformal field theories parametrized by the deformation parameter (coupling constant). 
The anomalous dimensions A depend on this parameter. In the model (1) this situation 
is realized when one introduces an Abelian disorder (see 0). In the opposite case, when 
/ 7^ 0, i.e. the OPE of O{z)O{0) contains O itself, there is a renormalization group (RG) 
flow of the coupling constant 



dlnA 

which means that the theory actually depends on the scale A. 

Let us now study the same problem in a case where the theory is deformed by the 
operator DD, which is truly marginal, because the OPE of D{z)D{0) does not contain the 
operator D itself (see Eg. (pOD . In this case we shall calculate the correlation function 



G{z;-f)= {A{z)B{0)exp{'y J d^xDD{x)) 

= E ^ / {Mz)B{0) DD{xi)....DD{xn))d'xi d^Xn (63) 

where A and B are some operators [Q and Q^, for example) with the correlation function 

G{z- 0) = (A(z)5(0)) 

Using the OPE 

1 2 I 1 2 1 2 I 1 2 

DD{x) A{y) = a " ^' A{y), DD{x) B{y) = & ^ ^' B{y) (64) 
\x y\ \x y\ 

one can find the first order in 7 correction to the correlation function which will be 



7 j {A{z)B{0) DD{x))d^x = (65) 

aj{A{z)B{0)) f d^x^-^P^^^ +bj{A{z)B{0)) f d\^-^P^ 
J \x — zY J 

where in both integrals we integrate over x between and z. Then it is easy to find the 

following logarithmic correction: 

Now one can consider the next order corrections and sum all of them using the same methods 
as in the case of conventional marginal operators (see fl^)- The result is 

G{z- 7) = G{z- 0) exp (^7^ In' kp) (66) 

which is different from the case of a conventional marginal operator when one has the first 
power of log in the exponent and not the third. The first power in the exponent introduces 
the power factor 

exp ((a + 6)7 In = |^|2(«+% (67) 

corresponding to the change in the anomalous dimension Ayi(A) = A^(0) — (0 + 6)7, and the 
behaviour of the deformed correlation function is still power-like. The are no logarithmic 
corrections after all. This is not true anymore with the logarithmic operator, when the 
correlation function cannot be written as a power at all. Thus we see that the correlation 
functions for operators which have non-trivial OPE with the logarithmic operator D (like 
our primary fields Q, for example) will have logarithmic corrections in the deformed theory 
- even in the absence of the RG fiow. 

hei A = Q and B = Q+ , then, according to Eq.(|3D, a = b = N'^. At 7 < the 
correlation function ( |66D decays faster then any power. At 7 > it increases faster than any 
power. In this case the approximation leading to Eq.(|66D breaks down when the correlation 
function begins to increase, that is at 

~ M-^exp[-^] (68) 
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We speculate that for 7 > the symmetry is broken and the finite density of states at e = 
is formed. This probably explains the finite DOS obtained numerically in disordered d-wave 
superconductors by Wheatley |J15 . 



VII. PROBABILITY DISTRIBUTION OF LOCAL DOS 



Now we shall calculate the distribution function of local densities of states. We can do 
it for a system of a finite size L. From Eq.(^ we know that in the case of zero frequency we 
have 



(p"(x)) = M"([Tr(Q + Q+)n 



-2A„ 



(69) 



The latter equality is valid in the leading order in A„ is given by Eq.(|l^) being the 
smallest conformal dimension in the operator product of n operators Tr(((5 + Q^))- For 
N > 2 A„'s are negative for n > 1. 

Let us imagine now that the result (|69D comes from a local distribution function of p: 



(N - Dri^ 

P(p)p"dp = A„exp{21nL[^^^ 



— ]} 



(70) 



where may contain powers of InL. The distribution function which reproduces this 
result is the famous log-normal distribution which is considered as a characteristic feature 
of disordered systems [^], |]T^, |T^: 

1 



P(p)=D(p)exp 



In L"^ 



(71) 
(72) 



where -D(p) is a smooth function of Inp which we cannot determine. 

From the fact that the frequency scales as 2,-2+2Ai ^^^^ conjecture that the distribution 
function of p{uj) is given by 

1 



P{puj) = D{p^) exp 



ln(l/cj^) 
3A^-2 



ln'(p.^~^ 
A{N-1 



2A^2 _ 1 



,7 



2m 



(73) 
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The authors of Ref. 0] have discussed this distribution without writing it down exphcitly. 

It is worth remarking that the log-normal distribution also appears in the Liouville theory 
(27). Indeed, according to Eg. (p2D , the conformal dimension of the operator exp{na(f)) is 
given by ^na{Q — na), i.e. it has the same quadratic n-dependence as A„ (|1^). Repeating 



the previous arguments we obtain for the vertex operator V = exp{a(f)) the log- normal 
distribution (0) with r] = 4a^, ( = aQ{3 — 2aQ). Since the logarithmic operators appear 
in this theory only when Q = 2a, we conclude that their presence is not directly related to 
multifractality of the target space. 

VIII. CONCLUSIONS 

In this paper we have demonstrated that in the general class of nonunitary critical models 
a new phenomenon takes place - the emergence of logarithmic operators associated with a 
special hidden continuous symmetry. The presence of this symmetry is intimately related to 
the fact that the order parameter of our model - the local density of states at e = - does 
not acquire a non-zero average. Since above mentioned features appear also in quantum 



gravity |10|, we anticipate a connection between quantum gravity and critical models with 



disorder |I8|. Our expectations are supported by the recently discovered similarities between 



the conventional localization theory and the Liouville theory [|l^, pO 



The physical meaning of the hidden symmetry in models with disorder as well as in 2d 
gravity, remains obscure. It is clear, however, that this symmetry should routinely appear 
in critical non-unitary theories where the Hamiltonian cannot be diagonalized (in our paper 
this fact is expressed in Eq.(^). 

It follows also from our work that, at least for the model in question, the replica approach 
is equivalent to the supersymmetric one. This is a pleasant fact. 

As we have said above, Eqs.(^) and ( |55D enable one in principle to reformulate the theory 
in terms of representations of the current algebra of the conserved current C{z). With this 
task being accomplished one can abanbon replicas and treat the theory axiomatically as it 
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is customary, for instance, in the theory of the standard WZNW model. At present this 
remains the biggest challenge. 
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Appendix 

In this Appendix we write down the relationship between conformal blocks in the replica 
and the supersymmetric (SUSY) representations. According to Ref. 0], the correlation 
functions in SUSY representation are products of correlation functions of the Gaussian 
model and the SUfc(N) WZNW theory with k = —2N. According to Ref. |0, the conformal 
blocks of the latter model are 

Wi°\x)= {l-x)-^F{l/N,-l/N,2;x) 

W^f (^)= -^^(1 - ^r'H^ + ^/N, 1 - 3; x) (74) 

The second solutions contain logarithms. The relationship between two representations is 
established by the identity 

x{l - l/N^)F{l + 1 - 3; x) + 2(1 - x)F{l/N, -l/N, 2; x) 

= {l-xN'^)F{l/N,-l/N,l-x) (75) 

using which one can write the expression for four-point function ([T8| ) in terms of conformal 
blocks either of replica or supersymmetric models. 
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